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The persistent current and charge stiffness of a one- 
dimensional Luttinger liquid on a ring threaded by a magnetic 
flux are calculated by Monte Carlo simulation. By chang- 
ing the random impurity potential strength and the electron- 
electron interaction, we see a crossover behavior between weak 
and strong impurity limits. For weak impurity potentials, in- 
teractions enhance impurity effects, that is, interactions de- 
crease the current and the stiffness. On the other hand, inter- 
actions tend to screen impurities when the impurity potential 
is strong. Temperature dependence of the persistent current 
and the charge stiffness shows a peak at a characteristic tem- 
perature, consistent with a recent single impurity study. 
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Recent technological progress novif allows us to produce 
nearly ideal low-dimensional systems, and it is becoming 
more important to understand how disorder and electron- 
electron interactions control physics in such systems. In 
one-dimensional systems, Giamarchi and Schultz |3] de- 
rived renormalization group (RG) equations for systems 
with weak but finite density of impurities, and showed 
that the critical impurity strength is zero for the Lut- 
tinger liquid parameter g > 2/3 and is finite for g < 2/3, 
where g > 1 and g < 1 correspond to repulsive and at- 
tractive interactions, respectively. Scalettar et al. [^ con- 
firmed a part of the phase diagram given by Giamarchi 
and Schultz, using quantum Monte Carlo (MC) simula- 
tion of the boson Hubbard model. Q] 

For finite-size systems, current can flow even in the 
presence of impurities for any g because of the finite size. 
One interesting example of a finite-size system is a one- 
dimensional ring threaded by magnetic flux, which cre- 
ates persistent current as first predicted by Byers and 
Yang and later by Biittiker et al. B. When an electron 
gas moves in a dirty ring, we expect the persistent cur- 
rent to be suppressed as the impurity potential becomes 
stronger B . What then will happen if electron-electron 
interactions are turned on 0? Krive et al. |g] calculated 
persistent current of Wigner crystal on a ring with a sin- 
gle impurity in certain limits. From their expressions 
we can see the current decreases as we increase particle- 
particle repulsive interactions. This study indicates that 
interactions effectively enhance the impurity strength, al- 
though their system has only one impurity. For a finite 
density of weak impurities, we can use RG equations ob- 
tained by Giamarchi and Schultz, placing a cutoff in the 
length scale since the system is finite. Then the effective 



impurity strength W is given by VF ~ L^~^/s^ where L is 
the system size and g is the Luttinger liquid parameter 
which becomes large as repulsive interactions increase. 
Here we can see again that interactions would enhance 
the effective impurity strength. 

There have been several exact diagonalization studies 
done for the ring problem |^. These calculations were 
performed in diffusive regime, i.e., for relatively weak 
impurities, and showed qualitatively the same results as 
the above analytical works. 

This interplay between interactions and impurities can 
be understood in the following way. When impurities are 
weak enough, the system still resides in the Luttinger 
liquid regime |In|. Luttinger liquids have charge-density 
wave (CDW) correlations which grow as the repulsive in- 
teraction becomes strong. Since a CDW would be easily 
pinned by impurities, stronger interactions make the sys- 
tem pinned more easily. Therefore the persistent current 
would decrease when interactions increase. 

Now the question is what will happen with strong im- 
purities. It seems clear that interactions would screen 
the impurities, which is opposite to what is expected in 
the weak impurity limit. When the impurities are very 
strong, particles are almost localized in the low poten- 
tial regions. If we turn on repulsive interactions between 
particles, it will move the particles away from each other, 
and therefore the interactions would effectively weaken 
the impurity strength. 

In this letter we will show the crossover behavior be- 
tween these weak and strong impurity limits. We have 
calculated the persistent current and charge stiffness of 
a spinless Luttinger liquid on a ring with a finite density 
of impurities, using Monte Carlo simulation. One ad- 
vantage of working with a Luttinger liquid instead of a 
microscopic lattice model is that we can separate out the 
lattice (periodic potential) effect, which is another impor- 
tant factor changing those quantities jll| , so we can con- 
centrate on the interplay between impurities and inter- 
actions only [Q. We also show the temperature depen- 
dence of persistent current and charge stiffness. These 
have peaks at certain temperatures, as previously argued 
for the single impurity problem [g|. 

The problem of Luttinger liquid on an impurity-free 
ring penetrated by a magnetic flux was solved by Loss 
]l3[ . He showed that the Lagrangian is given at T <C 
Trgx/2L by 
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id^ef + {N-i + 2<i)/M^dre, 



(1) 



where ^jl = t and x. The boundary condition of 9 field 
is 9(x + L,T + P) — 9{x, t) + 2n7r, where n is a winding 
number. The derivatives of 9 give the particle density 
and current: d^9 — —2ttj^. Defining 9 = 9 — 2nT7r//3, we 
obtain the partition funtion Z in the following form, 

Z = y /d07-l)(^-i)"cos(27rn^)e~*"'e-/^ (2) 



where C = {g^/ Sir) (8^9)'^, and Tq = 2/grLTr, In the fol- 
lowing calculation we assume the particle number N is 
odd, since we can get the physical quantities for the cases 
where N is even by using the relation Z{N even; (f>/4>o) = 
Z{N odd;cj)/(j}o + 1/2). Haldanc ||l|| calculated the re- 
lation of Qx and gr to the microscopic parameters of a 
spinless fermion system with nearest neighbor interac- 
tions, gr varies with density but does not change very 
much with interaction as long as the interaction is not 
too large, while gx changes with both density and inter- 
action. The sound velocity Vs is given by \/ gxl Qr- In the 
case of the noninteracting model, gx — ^/ gr — Vs = vp. 

The impurity term is written in the microscopic model 
as J dxViinpix)p{x). Following Rcf. [|5| we obtain 



/ Vi,np{x)pix) = / ^Vimpdx9 + J 



Vimp{x)poe 
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-h.c. + higher orders, 



(3) 



where po is average density. The first term is forward 
scattering and the second represents backward scattering. 
Using 9 instead of 9 we obtain 



C^^{d^,9f + ^V,,^p{x)dx9 

OTT ZTT 

-t-2/9oV^mp(a;) cos(6' -I- 27rpoa; + 2n7rT//3), 



(4) 



where we have omitted the tilde from 9 and we used 
kp = TTpo- Since the forward scattering is irrelevant pj, 
we took into account only the first and third terms in our 
calculations. 

In order to do MC simulations, we do not want to have 
the Aharonov-Bohm phase factor cos(27r(/)/(/)o) in Eq.(2) 
because it can give a negative weight. Fortunately we can 
express persistent current for finite in terms of (n'^) 1^=0 
as will be shown below, so that we just need to perform 
MC simulations at (/> = which causes no negative weight 
problem. 

Here we show how to calculate persistent current at 
finite (f> using quantities ai (f> — 0. The persistent current 
/ on a ring is given by I{(p) = —{l//3)d\nZ{(j))/d(p = 
— (27ri//3(/)o)(n). We expand I{ip) in a Fourier series: 
H4>) = Z]„-f"Sin(27rn0/(/>o), in other words, (n)(0) = 
Pi{(j)o/2n) ^„ /„ sin(27rn0/(/>o). Taking derivatives of (n) 
with respect to and setting = 0, we obtain. 



d^'+Hn), , ,,,,Ao,,,- 
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^u^o = (-f/?^(£rE^" 



,2/c+l 



(5) 



and d {n)/d4> \^^Q — On the other hand, from the 
definition of (n) we have 



0(j) 00 



(6) 



^U=o = (^)^((n^>-3(n2)^)U=o, (7) 



and d {n) / d4> \^=Q = 0. Consequently, we have 

27r 
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(8) 
(9) 



If we assume /„ for n > 2 are negligibly small, we can 
solve Eq. (||) and (^ for h and I2, 

/i = ^(-(O+3(n2)2+4(n2))|^=0, (10) 

h = ^{-{n'')+HnY + {n^))U=,. (11) 

Therefore we can calculate approximately the first and 
the second harmonics of persistent current from expec- 
tation values (n^'^) at = 0. The charge stiffness 
is proportional to the mean-square winding numbers: 

Ps = Ld^F/dq^\^o = {L/P){n^)\^=o- 

Since our Lagrangian is invariant under the global 
transformation 9—^9 + 27rfc with k being an integer, 
the 0's could go far away from the initial values during 
the updating processes of MC simulations. In order to 
prevent this behavior which is inconvenient to the simula- 
tions, we added a mass term (m/2)9'^ to the Lagrangian. 
After some experiments we found m — 0.02 is appro- 
priate and we used this value in the following numerical 
calculations. We chose gr = l/-\/2 which corresponds ap- 
proximately to quarter filling in the microscopic model of 
spinless fermions with nearest neighbor interactions [[l4| . 
We discretized the system in spatial and imaginary time 
directions, so that the system size and the inverse tem- 
perature are characterized by integers L and L,-, respec- 
tively. In the simulations we fixed the system size L = 8 
and the particle density po = 0.25. We took 10,000 warm- 
up steps and 60,000 MC steps for measurements. The im- 
purity average was taken over 300-600 realizations using 
the flat probability bounded by W: —W < V\mp(x) < W. 
Since there is a term (n^)^ in Eq. (|lO| ) and ( [ll[ ) we made 
two independent MC runs for each impurity realization 
in order to minimize the statistical errors in the impurity 
average. 

Figures |l| show the first harmonics /i and charge stiff- 
ness Ps as functions of impurity potential strength W 
for different interaction strengths. Temperature is fixed 
to 1/I/T = 1/24. We can see the crossover behavior be- 
tween weak and strong impurity limits. As known previ- 
ously, persistent current and charge stiffness are indepen- 
dent of interactions when there is no impurity {W = 0) 



P, pl]Jl^ . When W is small, both Ii and ps are sup- 
pressed as gx increases, i.e., as repulsive interactions in- 
crease, whereas they increase with g^ in the strong im- 
purity region {W > 2). We mentioned earlier that the 
crossover should occur because the system resides in the 
Luttingcr liquid regime for weak impurities whereas the 
particles are localized with strong impurities. There- 
fore we believe the crossover occurs when the localization 
length reaches the system size. 

The temperature dependence of Ii and ps are shown 
in Figs. 2 and 3. We can see that Ji and Ps have peaks 
at certain temperatures. This peak behavior comes from 
two competing effects: depinning by thermal activation 
and temperature destruction of quantum coherence. The 
former effect is important at rather low temperature but 
the latter dominates at higher temperature. This be- 
havior has been seen in a single impurity problem for a 
Wigner crystal M, and here we find the essentially the 
same feature in the Luttingcr liquid in the presence of 
finite density of impurities. The peak position T* is 
expected to be determined by those competing effects, 
that is, T* should be controlled by the effective impu- 
rity strength. As discussed above, the effective impurity 
strength in the weak impurity region is enhanced by in- 
teractions. Therefore T* should be larger for larger g^. 
Actually this is what we can see in Figs, g (a) and (b). 
Figures. S (a) and (b) show that, in strong impurity 
region, T* is shifted to the lower side as gx increases, be- 
cause interactions would reduce the impurity strength. 

Another interesting feature is the role of interactions 
at rather high temperature. For weak impurities, inter- 
actions always suppress the persistent current and the 
stiffness in the whole temperature range. When we have 
strong impurities, however, interactions would screen the 
impurities only at low temperature, and would suppress 
the current and the stiffness at high temperature. This 
could be because electrons are already delocalized by the 
thermal activation at high temperature but the CDW 
correlation still remains, in other words, the localiza- 
tion length becomes shorter than the system size above 
some temperature although the CDW correlation length 
is large. Therefore the role of interactions in strong im- 
purity region at high temperature is expected to be the 
same as in weak impurity region. 

In summary, we have reported on the first MC simu- 
lation of the ID ring problem to see the interplay be- 
tween disorder and electron-electron interactions. We 
have shown a crossover behavior from the weak-impurity 
to strong-impurity region. When the impurities are weak, 
interactions try to suppress the persistent current and 
charge stiffness. On the other hand, interactions enhance 
the current and the stiffness in the presence of strong 
impurities. In other words, interactions increase the ef- 
fective strength of weak impurities, and decrease the ef- 
fective strength of strong impurities. We gave a physical 
picture of these opposite behaviors, and we believe the 
crossover occurs when the localization length becomes or- 
der of the system size. The temperature dependences of 



persistent current and charge stiffness are also calculated. 
We showed they have peaks at certain temperatures and 
the peaks shift as we change interaction strength. The 
shifts are consistent with the change of the effective im- 
purity strength by interactions. At high temperature the 
thermal activation would screen the impurity effect and 
interactions act in the same way as they do with weak 
impurities. 
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FIG. 1. (a) Impurity and interaction dependence of the 
first harmonics h of persistent current, g^ is a renormalized 
parameter which is an increasing function of the microscopic 
particle-particle interactions, (b) The same dependence of 
charge stiffness ps- 



FIG. 2. Temperature dependence of 7i (a) and ps (b) in 
the weak-impurity region (W=l). The dashed lines are for a 
Luttinger liquid with no impurity. 



FIG. 3. Temperature dependence of Ji (a) and ps (b) in 
the strong- impurity region (W=3). 
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